Phonon-Magnon coupling in CoF 2 investigated by time-of-flight neutron spectroscopy 
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We report the results of inelastic neutron scattering investigation on the model antiferromagnet 
C0F2 by time-of-flight neutron spectroscopy. We measured the details of the scattering function 
5(Q, uj) as a function of temperature with two different incident neutron wavelengths. The tempera- 
ture and Q dependence of the measured scattering function suggests the presence of magnon-phonon 
coupling in almost all branches. The present results are in agreement with the strong magnetoelastic 
effects observed previously. 



PACS numbers: 75.25,+z 

The lattice dynamics of a crystal is normally assumed 
to be independent of the spin dynamics of the magnetic 
moments associated with the atoms and ions of the crys- 
tal. Similarly while describing the motions of the spins in 
a magnetically ordered systems it is usual to assume that 
the equilibrium positions of the magnetic ions are frozen. 
The approximation gives a fairly well description of the 
magnetic system and is in the same spirit as the adiabatic 
or Born-Oppenheimer approximation which is normally 
used for band structure calculations. In this approxima- 
tion the phonon and magnons are separate entities with 
their own dispersion relations and the Hamiltonian of the 
system separates out into a lattice and a magnetic part. 
In reality however there always exists a coupling term 
that involves displacements both in the positions of the 
atoms in the crystal and in the orientations of the spin 
vectors associated with these atoms. Thus the collective 
excitations are neither pure lattice waves nor spin waves 
but rather magnetoelastic waves. Such magnetoelastic 
waves are predicted to exist theoretically long time agcffl 
With the advent of enormous progress in first-principle 
calculations it is now possible to do phonon calculations 
of magnetic solids in the magnetically ordered phase by 
taking spin degrees of freedom into account. The dif- 
ference between the results of the calculations done by 
taking and not taking magnetic degrees of freedom into 
account can be often very significant as has been amply 
demonstrated by such recent calculations on Fe-arsenide 
compounds^. It is possible to study magnetoelastic ef- 
fects experimentally by using neutron scattering tech- 
niques. Our recent neutron diffraction investigations^^ 
have shown strong static magnetoelastic coupling in tran- 
sition metal difluorides MnF2, FeF2, C0F2 and MF2. 
There exist indirect evidences for dynamic magnetoe- 
lastic or phonon-magnon coupling in FeF2, C0F2 and 
MF2 56 . However, a systematic study of the phonon- 
magnon coupling in this relatively simple system is still 
missing. 

C0F2 belongs to the family of transition-metal di- 
fluoride, which has been the subject of intensive 
investigations^^. C0F2, along with other transition- 
metal oxides MnF2, FeF2 and N1F2, crystallize with the 




FIG. 1: The antiferromagnetic structure adopted by 3d tran- 
sition metal diflurides viz. C0F2 with the rutile type crystal 
structure. The Co and F ions are labeled. The arrows on the 
Co represent moment direction of the transition metal ions 
below the Neel temperature 31 . 



tetragonal rutile- type structure in the P^/mnm space 
group. However the magnetic properties of C0F2 are 
more complex than those of isomorphous MnF2, because 
the Co ion has unpaired angular momentum that plays 
an important role in determining its magnetic proper- 
ties. CoF 2 orderPJ^Sl below T N w 39 K with an an- 
tiferromagnetic structure^, shown in Fig. [I] with the 
propagation vector k = 0. The magnetic moments of Co 
ions at the corner (000) positions of the tetragonal unit 
cell are all parallel to the c-axis whereas those of the Co 
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FIG. 2: (Color online) Temperature evolution of the scat- 
tering function S(Q,u) measured from CoF2 using incident 
neutron wavelength of A = 1.5 A. 



ions at the (|||) positions are oppositely oriented. The 
presence of an orbital moment in C0F2 makes it particu- 
larly interesting for the study of magnetoelastic coupling 
in this compound. We recently investigated^ the hy- 
perfine interaction in C0F2 by high resolution neutron 
spectroscopy and concluded that the presence of the un- 
quenched orbital moment of the Co ion in C0F2 leads to 
its anomalous behavior compared to that of other Co- 
based compounds. 

We have performed inelastic neutron scattering inves- 
tigations on C0F2 using the thermal time-of-flight neu- 
tron spectrometer IN4C at the Institute Laue-Langevin 
in Grenoble. About 5 g of C0F2 powder sample was put 
inside a cylindrical Al sample holder that was fixed to the 
cold tip of the sample stick of a standard orange cryostat. 
In order to cover an extended Q-range, relevant to the 
present study, and to gain in energy resolution two inci- 
dent wavelengths of A^ = 1.1 and 1.5 A, respectively, were 
used. This allowed the Stokes spectrum to be measured 
at low temperature over a broader dynamical range. The 
data analysis was done using ILL software tools 
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FIG. 3: (Color online) Temperature evolution of the scatter- 
ing function S(Q, uj) from CoF2 measured with the incident 
neutron wavelength of A = 1.1 A. 



The present investigations have been performed with 
unpolarized neutrons on powder samples of C0F2. The 
present method is therefore somewhat limited because it 
gives only the powder average results and does not give 
information of detailed directional Q dependence. Also 
the unpolarised neutron scattering cannot distinguish di- 
rectly whether the scattering is of magnetic or structural 
origin. However and since C0F2 contains magnetic Co 
ions, we can therefore expect inelastic magnetic scat- 
tering due to the interaction of the neutron spin with 
the spin and orbital magnetic magnetization of the C0F2 
both below and above the magnetic ordering temperature 
T/v ~ 39 K. Also the neutron is sensitive to the struc- 
tural excitations or phonons which exist at all tempera- 
tures. Thus this technique will measure both excitations 
and even highlights their hybrid feature described before. 
Magnetic and structural excitations have fortunately dif- 
ferent temperature and Q dependence. Strong and sharp 
inelastic magnetic excitations are expected only below 
the magnetic ordering temperatures. At higher temper- 
ature, the magnetic correlations give broad quasielastic 
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FIG. 4: Results of the fit of the low-Q data (integrated over 
Q in the range from 0.9 to 3.56 A -1 ) from C0F2 in the low- 
energy range (A* = 1.5 A). Two Gaussian peaks have been fit- 
ted to the elastic peak at zero energy and the inelastic peak 
which is at about E — 8 meV at T — 2 K. The energy of 
inelastic peak decreases with increasing temperature and be- 
comes zero at Tn = 38 K giving rise to quasielastic scattering 
which has been fitted with a Lorentzian function. 

scattering instead. The intensity of the magnetic excita- 
tion peaks decrease with the momentum transfer Q. The 
scattering due to structural excitations or phonons have 
different temperature and Q dependences. The phonon 
intensity in general increases with increase of the momen- 
tum transfer Q. Of course intensities of magnons and 
phonons are also governed by the magnon and phonon 
structure factors and also polarisations which depends 
on the details of the magnetic and crystalline structures. 
However the differences in temperature and Q depen- 
dence for magnons and phonons can be conveniently used 
for the sake of identification in order to separate them. 

The cross section for neutron scattering by magnetic 
systems is discus sed in details in several text books on 
neutron scatterin g 33 * 34 *. The expression for neutron scat- 
tering cross sections are very complex when the ions have 
both spin and orbital angular momentum. It can be 



FIG. 5: Results of the fit of the low-Q data (integrated over Q 
in the range from 0.9 to 3.56 A -1 ) from C0F2 in the higher- 
energy range (A* = 1.5 A). At low temperatures below Tn = 
38 K, three peaks could be identified and they were fitted 
with three Gaussian functions. At higher temperature only 
a broad peak could be identified and has feen fitted with a 
single Gaussian function. 



greatly simplified if the momentum transfer is less than 
the reciprocal of the radius of the magnetic shell and is 
given by 

d 2 a _ k f 1 /1.91e 2 \ 2 
dQduj ~ Yi2^ V 2mc<2 ) 

/CO 
<K(t,0).K(7,f)> 

x exp [iQ.(R(i) - R(j))] exp(—iui)dt, 

where /(Q) is the form factor of the ions, K(j, t) is the 
part of the magnetic moment operator, L + 2S, which is 
perpendicular to Q, for the ion i ar R(i) at the time t, Iq 
and kj are the wavevectors of the incident and scattered 
neutrons. 

The neutron scattering cross section of a phonon mode 
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FIG. 6: (Color online) (a) Temperature variation of the en- 
ergy and the half-width at half-maximum (HWHM) of the 
low-energy inelastic magnetic peak of C0F2. (b) Temperature 
variation of the energy of the high-energy inelastic magnetic 
peak of CoF 2 (A; = 1.5 A). 



(qj) is given by 



d 2 a 



k f 1 



dftduj ki uj(qj) 



b k Q.e(k,qj) 



(2) 



2M k 



x exp[z(Q - q).R(fc)] eM-W k (Q))\ 2 

n(qj) s 



n(qj) + Y 



where exp(— W k (Q)) is the Debye- Waller factor for the 
fcth atom, n (qj) is the population of the mode and is 
taken for neutron energy loss, b k is the scattering length 
of the kth atom in the unit cell, M k its mass, and e(fc, qj) 
is the eigen vector in the normal mode (qj). Note that 
everything we stated in the previous section about the 
different Q dependence of the magnons and phonons, fol- 
lows directly from equations ([I]) and ([2| 



FIG. 7: Q dependence of the intensity for the inelastic peaks 
at E = 7.3, 14.0, 23.5, 31.6 and 43.0 meV at T = 2, 20, 40, 
70 and 100 K (A* = 1.1 A). 



Both magnetic and structural excitations have been 
investigate d 27 * 28 * by inelastic neutron scattering on C0F2 
single crystals. The magnetic properties of C0F2 are 
much more complex than those of MnF2 having the same 
rutile crystal structure and also the same antiferromag- 
netic structure. Crystal field calculations of Gladney^H 
shows that the L = 3,S = 3/2 atomic ground state is 
split by the octahedral crystal field to give an orbital T4 
triplet as the ground state whose properties can be de- 
scribed by an effective I = 1 operator. The rhombic and 
tetragonal crystal field distortions and the spin-orbit cou- 
pling causes this state to split into six Kramers doublets. 
The molecular field in the antiferromagnetic phase causes 
further splitting. Labeling the ground state as A and 
first three excited states as B,C and D three crystal-field 
excitations have been labeled 27 as A-B, A-C and A-D 
transitions. Alternatively these states are also labelecP 
as 0, 1, 2, 3. The corresponding excitations are labeled 
as - 1, - 2 and 0-3. The dispersions of these exci- 
tations have been measured from single crystal samples 
along [100] and [001] directions. However, the anisotropy 
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in dispersion of these excitations is found to be rather 
small. The strongest excitation is A-B which is at about 
7.4 meV for q — 0. The other two excitations are A-C 
and A-D that both lie at about 24-25 meV. 

Fig. [2] shows the temperature evolution of the scat- 
tering function S(Q,uj) measured with incident neutron 
wavelength of A = 1.5 A and Fig. [3] shows temperature 
evolution of the scattering function S(Q,u) measured 
with A = 1.1 A. At T — 2 K we observe two inelas- 
tic peaks at about E = 8 meV and E — 24 meV. The 
energy of these peaks decreases with increasing temper- 
ature. However whereas the energy of the first peak at 
E = 8 meV becomes almost zero at T/v = 38 K, the 
energy of the second peak decreases and then attains a 
finite value at T/v and stays constant at about 19 meV 
at higher temperatures. Fig. [4] shows the results of the 
fit of the low-Q data from C0F2 in the low-energy range 
in the left two panels. Two Gaussian peaks have been 
fitted to account for the elastic peak at zero energy and 
the inelastic peak which is at about E = 8 meV at T = 2 
K. The energy of the inelastic peak decreases with in- 
creasing temperature and becomes zero at T/v = 38 K 
giving rise to quasielastic scattering, which has been fit- 
ted with a Lorentzian function. The results of the fit 
with Gaussian functions to the higher-energy peaks are 
shown in Fig. [5] At low temperatures below T/v = 38 K, 
three peaks could be identified and they were fitted with 
three Gaussian functions. At higher temperatures only a 
broad peak could be identified and has feen fitted with a 
single Gaussian function. Fig. [6] (a) shows the temper- 
ature variation of the energy and the half-width at half- 
maximum (HWHM) of the low-energy inelastic magnetic 
peak of C0F2 . In Fig. [6] (b) the temperature variation of 
the energy of the higher-energy inelastic magnetic peak 
of C0F2 has been shown. The temperature variation of 
the low energy inelastic peak which is at about 8 meV 
at 2 K shows clearly that this peak is due to magnetic 
excitations. The energy of this peak decreases continu- 
ously and becomes almost zero at T/v = 38 K giving rise 
to quasielastic scattering of Lorentzian form. The higher 
energy peak at about 24 meV on the other hand seems 
to be a hybrid peak due to the magnon-phonon coupling. 
The peak persists at temperature higher then T/v and its 
energy remains constant at about 19 meV. 

Fig. [7] shows the Q dependence of the intensity for 
the inelastic peaks at E — 7.3, 14.0, 23.5, 31.6 and 43.0 
meV at T = 2, 20, 40, 70 and 100 K. Fig. [7] (a) shows 
the Q dependence of the peak at E = 7.3 meV for sev- 
eral temperatures. The intensity of the peak at E = 7.3 
meV at T = 2 K decreases continuously with increasing 
Q and then becomes more or less constant at higher Q. 
This indicates that the inelastic peak is of mostly mag- 
netic origin. This is also supported by its intensity and 
Q dependence at T = 100 K. However even at T = 100 K 
the intensity is not zero. This suggests that there exists 
some phonon contribution to this peak as well. The Q 
dependence at T = 100 K shows that the intensity in- 
creases slightly at higher Q and this also confirms that 



there exist some phonon contribution to this peak. Fig. 
[7] (b) shows the Q dependence of the inelastic peak at E 
= 14 meV at several temperatures. The intensity of this 
peak first decreases with increasing Q and then increases 
very much at higher Q. This suggests that the peak has 
a hybrid phonon and magnon character. The Q depen- 
dence of such a hybrid peak is also expected to be more 
exotic. Fig. [7] (c) shows the Q dependence of the inelastic 
peak at E = 23.5 meV. Its Q dependence is very similar 
to the peak at E = 7.3 meV shown in Fig. [7] (a). Even 
its temperature dependence is very similar. The inten- 
sity of the peak at all Q is higher at lower temperature. 
The decrease of its intensity with increasing Q suggests 
that it is magnon like. However like the peak at E = 7.3 
meV the peak at E = 23.5 also shows slight tendency of 
increasing in intensity at higher Q suggesting that the 
peak has some phonon like character. Fig. [7] (d) shows 
the Q dependence of the inelastic peak at E = 31.6 meV. 
It looks very similar to that at E = 14 meV and con- 
sists clearly of both phonon and magnon contributions. 
Lastly the Q dependence of the inelastic peak shown in 
Fig. [7] (e) shows that this peak is almost pure phonon 
like because its intensity increases with Q and tempera- 
ture. We have fitted the Q dependence of the intensity 
for the peaks at 14 and 31.6 meV assuming them to con- 
sist of both magnon and phonon contributions and also 
assuming that the Q dependence has oscillatory cosine 
wave form due to the periodicity of the lattice. We fitted 
the intensity by the equation 

I(Q) = ai + a 2 Q + a 3 Q 2 + a 4 cos \—Q + a 6 J (3) 

W J 

which consists of a polynomial plus an oscillatory cosine 
function. Here ai,<22, and as are coefficients of a sec- 
ond order polynomial and is the amplitude and a$ is 
the periodicity of the oscillatory part superimposed on 
the polynomial function and is the phase. This is a 
much simplified function but since it contains term pro- 
portional to Q describing Q dependence of the magnetic 
intensity and a term proportional to Q 2 describing the 
Q dependence of the phonon intensity and has also a su- 
perimposed oscillatory term describing the effect of lat- 
tice periodicity, it describes the intensity variation rather 
well. The fitted periodicity is in agreement with the aver- 
age lattice periodicity. Since the fits for different temper- 
ature in Fig. [7] (b) is not clearly visible we show a typical 
fit at T = 2 K separately in Fig. [7] (f). The fit looks 
rather convincing indeed. In fact the Q dependences of 
the intensities for all five energies shown in Fig. [7] (a-f) 
show oscillations but the statistics and also the resolu- 
tion of the data shown in Fig. [7] (a), [7] (c) and [7] (e) are 
not good enough for reasonable fits. 

In conclusion we investigated the phonon and spin dy- 
namics of C0F2 by inelastic neutron scattering on powder 
samples by time- of- flight neutron spectroscopy with un- 
polarized neutrons. The temperature and the Q depen- 
dence of the intensity of the inelastic signals enabled us to 
identify their phonon, magnon or hybrid origin at least to 
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a first approximation. We note that there exist consider- 
able phonon-magnon interaction in this simple antiferro- 
magnetic system also evidenced by the strong magnetoe- 
lastic coupling observed by neutron diffraction 3 . How- 



ever, in order to be certain about the hybrid nature of 
these inelastic peaks and evaluate quantitatively phonon- 
magnon interaction strength, polarized neutron scatter- 
ing investigation is desirable. 
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